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ABSTRACT 

The null hypothesis that the three dimensional power spectrum measured from the 
APM Survey is consistent with the one dimensional power spectrum measured from 
the pencil beams surveys of Broadhurst et al. and of Szalay et al. is tested. The 
external estimates of the mean power that we make are sensitive to details of the 
model for the survey geometry and to the assumed level of the distortion of the 
pattern of galaxy clustering caused by peculiar motions. We hnd that the measured 
3D clustering of galaxies can account for the presence of peaks in the one dimensional 
power spectrum, but is less successful in recovering the detailed appearence of the 
observations. We find no strong evidence for any additional large scale structure in 
the deep pencil beams beyond that recovered from the APM Survey. We conclude 
therefore that it is unlikey that large scale structure can be responsible for the steep 
local number counts of bright galaxies. 

Key words: surveys-galaxies: clustering -dark matter - large-scale structure of Uni- 
verse 



1 INTRODUCTION 

1 The claim by Broadhurst et al. (1990, hereafter BEKS) that 
the universe contains significant structure on scales larger 
than ~ 100h~ Mpc has provoked much controversy. Based 
upon the analysis of several deep, narrow angle or pencil 
beam redshift surveys of galaxies in the directions of the 
Galactic Poles, BEKS noted a striking periodicity in the 
pair counts of the galaxies. This periodicity is revealed as a 
spike in the one dimensional power spectrum of the radial 
galaxy distribution at k ~ 0.05/iMpc -1 , corresponding to 
a wavelength of A ~ 128/i _1 Mpc. Using estimates of the 
mean power from the observed power spectrum (internal) 
and using a model for galaxy clustering combined with a 
window function for the survey (external), BEKS attatched 
a high significance to the peak. The adopted null hypothesis, 
namely that the two point correlation function of galaxies 
on small scales has a power law form f(r) = (r/ro) 7 (Davis 
& Peebles 1983), with a truncation at r = 30/i _1 Mpc and 
random correlations on larger scales, was ruled out. 

In a comprehensive paper, Kaiser & Peacock (1991 - 
KP) argued that the mean power had been underestimated, 
which would reduce the significance of the power spectrum 
peak. Using a simple Poisson clumps model, KP demon- 
strated that the power spectrum at higher wavenumbers 
could be suppressed by peculiar velocities and redshift er- 
rors, leading to a low internal estimate of the mean power. 
Furthermore, with more detailed modelling of the BEKS 
survey geometry and radial selection, KP showed that the 
external estimate of the power could be a factor of two larger 



than the figure obtained by BEKS. Note that KP employ 
a different null hypothesis, in that the observed small scale 
clustering of galaxies is extended to all scales using the same 
power law. 

Szalay et al. (1993) remark that assessments of the sig- 
nificance of the peak in the BEKS power spectrum based 
upon external estimates of the mean power are unreliable. 
Such estimates are sensitive to assumptions made about the 
clustering of galaxies on large scales and about the survey ge- 
ometry. This point is re-enforced by Luo & Vishniac (1993), 
who using the same model for galaxy clustering but a slightly 
different window geometry and selection function to KP, ob- 
tain a value for the mean power that is intermediate between 
that of Szalay et al. (1993) and KP. 

The original BEKS data has now been supplemented 
by additional pencil beams giving a wider angular coverage 
(Szalay et al. 1993, Broadhurst et al. 1995). These authors 
claim that the aliasing of small scale power to large scales is 
reduced for the extended survey geometry, which has a nar- 
rower window function in fc-space. In addition, there are new 
claims for excess power on scales of A ~ 100/i _1 Mpc, com- 
pared with standard models of structure formation such as 
the Cold Dark Matter scenario and its variants, from power 
spectrum analysis of the slices of the Las Campanas Red- 
shift Survey (Landy et al. 1996), though the longest baseline 
of this survey is shorter than that of the BEKS data. 

Measurements of the three dimensional power spectrum 
on scales A > 100/i _1 Mpc are now available from large red- 
shift surveys (e.g. Feldman, Kaiser & Peacock 1994, Tadros 
& Efstathiou 1996) and from the deprojection of the clus- 
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tering information in angular catalogues such as the APM 
Survey (Peacock 1991, Baugh & Efstathiou 1993, 1994). 

In this paper we shall ask the question is the three di- 
mensional power spectrum measured in real space from the 
APM Survey (Baugh & Efstathiou 1993, 1994) consistent 
with the one dimensional power spectrum of the BEKS sur- 
vey and the new data reported by Szalay et al. (1993) and 
Broadhurst et al. (1995). The models adopted for the sur- 
vey geometry of the observations are outlined in Section ^| 
with the effects of including redshift space distortions dis- 
cussed in Section []. The significance of the observed peaks 
are assessed by comparison with our external estimates of 
the mean power in Section [|. 



2 PENCIL BEAM WINDOW FUNCTIONS 

The projected one dimensional power spectrum is a con- 
volution of the three dimensional power spectrum with the 
Fourier transform of the window function of the survey: 



Pio(k) 
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(1) 



where I(k, y) depends upon the particular model adopted 
for the radial and angular selection of the survey. We shall 
discuss three survey geometries: 

A: Uniform cylinder: a disk of fixed comoving radius R, 
neglecting the radial selection function, with cylinder length 
L given by the highest redshift galaxies observed. Following 
BEKS, we use a cylinder with radius R = 3/i _1 Mpc and 
length L = 2034/i" 1 Mpc. 

B: A conical window, with fixed opening angle and a model 
for the survey selection function. In order to reproduce the 
deep NGP and SGP pencil beams of BEKS, we use a conical 
beam with opening angle 10 arcminutes and a magnitude 
limit of bj < 22.5. 

C: Multiple conical windows distributed at random within 
some larger solid angle. We choose 21 beams with opening 
angle 15 arcminutes distributed at random within a cone 
of opening angle 5°, with galaxies in the magnitude range 
17 < bj < 20.5, as described by Szalay et al. (1993). 

For a uniform cylinder, to second order in k the win- 
dow function takes the form (following BEKS, KP equation 
(3.10)) 



I(k,y) = T exp(- 
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To calculate the window function for a conical survey 
including the effects of a radial selection function, we follow 
the derivation given in Appendix Al of KP, substituting the 
redshift distribution of galaxies per steradian dj\f /dz used 
by Baugh & Efstathiou (1993): 

2zi 



cUV /dz = exp(- 



(z/zc) 1 - 5 ), 



(3) 



where the median redshift is given by 
z m (b.j) = 1.412z c 



= 0.016(6,/ 



17) 15 +0.046 



(4) 
(5) 



This parametric form for the redshift distribution was cho- 
sen to fit both the Stromlo-APM survey (Loveday et al. 
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Figure 1. The Fourier transform of the window functions for the 
survey geometries considered here. The solid line shows a uniform 
cylinder, the dashed line shows a conical window with a radial 
selection function and the dotted line shows the window function 
for multiple pencil beams. 



1992) and the fainter surveys of Broadhurst et al. (1988) 
and Colless et al. (1990, 1993) (c/ Figure 1 of Baugh & Ef- 
stathiou 1993). 

With this model for the selection function, the window 
function is given by: 



i%y) 
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(6) 



where the angular window function is normalised to unity. 
For a single beam, the Fourier transform of the angular se- 
lection function is given by 



W{kff) = exp(-(fc0) 2 /8), 



(7) 



whilst for N beams each of opening angle 8\ placed at ran- 
dom within a solid angle of opening angle 82, the angular- 
selection is given by (equation 4.5 KP): 

W(k) — W(k8i)(W(k9'2) + l/N), (8) 

with W(k8) given by equation 

The form of the window function Fourier transforms 
given above are plotted in Figure This Figure shows the 
large effect that the adopted survey geometry can have upon 
the deduced mean power; increasing the number of pencil 
beams damps the aliasing of power from small scales. 



3 REDSHIFT SPACE DISTORTIONS 

The pattern of galaxy clustering measured in redshift space 
is distorted by the peculiar motions of galaxies. Coherent 
bulk flows lead to a boost in the measured power in three 
dimensions on large scales (Kaiser 1987), whilst the viri- 
alised motions of galaxies in groups and clusters results in 
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Table 1. Mean power for survey geometries in real space 
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Figure 2. The Fourier transform of the window functions for the 
Case C, multiple beams. The solid line shows the window function 
in real space. The dotted line shows the window function when 
redshift distortions characterised by f3 = 1.0 and a v = lOOOkms - 1 
are included; the dashed line shows the window function for /3 = 
0.5, a v = SOOkms" 1 



a reduction in the power on small scales. Hence, in order to 
account properly for the effects of peculiar velocities when 
calculating the projected ID power, it is necessary to model 
both the small and large scale phenomena. 

A simple approach is to assume that galaxy velocities 
are given by linear theory on large scales plus an uncorre- 
cted, random velocity dispersion which dominates on small 
scales (e.g. Peacock & Dodds 1994; Cole, Fisher & Weinberg 
1995). The redshift space power spectrum Ps(fc,p) is then 
given in terms of the real space power spectrum Pfl(fc) by 



Ps(k,n) = Pj?(fc)(l + /3/x 2 ) 2 /(l + (ka vf i/H ) 2 /2) 2 



(9) 



where (3 = tt 6 /b, Q is the density of the universe in units of 
the critical density, b is the bias between fluctuations in the 
galaxy distribution and the underlying mass distribution, er„ 
is the ID velocity dispersion and /i = k ■ z/|k| is the cosine 
of the angle between the line of sight and the wavevector. 
We have assumed that the small scale velocities follow an 
exponential distribution as indicated by N-body simulations 
(Park et al 1994). 

Extending the results of Section 4.3 of KP the projected 
power including the effects of redshift space distortions on 
small and large scales is given by: 



Pio(fc) 
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The effects of including redshift space distortions on 
the shape of the window function in Case C, multiple pencil 
beams, is shown in Figure |j| Estimates of the parameter 
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Mean power 
±1(7 
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Ppeakl <P> 
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0.014 ±0.002 


15.2 


(B) 


bj < 21.5 
bj < 22.5 


0.065 ±0.009 
0.047 ±0.007 


3.3 
4.5 


(C) 


a = 
a = 2 


0.0085 ±0.0016 
0.0106 ± 0.0020 


10.0 
8.0 



Table 2. Mean power for geometry C in redshift space 
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Ppeak / ^ F > 


0.5 


500 


0.0107 ±0.0021 


7.9 


0.5 


1000 


0.0094 ± 0.0018 


9.1 


1.0 


500 


0.0127 ± 0.0024 


6.7 


1.0 


1000 


0.0150 ± 0.0028 


5.7 



(3 show a large spread in values. Determinations using the 
Stromlo/APM and APM Surveys by Loveday et al. (1995) 
and Baugh (1996) find ft ~ 0.5; a careful analysis of the 
errors by Tadros & Efstathiou (1996) shows however that 
values as large as /3 ~ 1 cannot be ruled out by the present 
generation of redshift surveys. Tadros & Efstathiou find one 
dimensional velocity dispersions in the range a v — 300 — 
700kms _1 in N-body simulations of the standard Cold Dark 
Matter model and its variants. 



4 PEAK SIGNIFICANCE 

We calculate the mean ID power spectrum at k = 
0.049/iMpc" 1 for each survey geometry, using the real space 
three dimensional power spectrum recovered in each of four 
zones into which the APM Survey was split by Baugh & 
Efstathiou (1993), with the results given in Table [l]. For ge- 
ometry B, we compute the mean expected power for two 
apparent magnitude limits. In case C for multiple pencil 
beams, we use two estimates of the real space power spec- 
trum in 3D. The parameter a describes the evolution of 
clustering with redshift; for a — clustering is fixed in co- 
moving co-ordinates, for a — 2 clustering evolves according 
to linear perturbation theory. 

In the case of a uniform cylinder, we note that the APM 
power spectrum gives a mean power that is lower than that 
obtained with a £(r) = (r/5) 1 ' 8 power law extended to all 
scales by 30% (KP find (P) = 0.019). Furthermore, if we 
truncate the APM power spectrum at A = 30/i _1 Mpc, the 
mean power for a cylinder falls by a factor of two. More 
realistic modelling of the survey window function boosts the 
mean power by a factor of 3 — 5, confirming the claim of KP 
that there is little evidence for periodicity in the deep NGP 
+ SGP beams alone. 

The final column in Table [l] gives the ratio of the ob- 
served peak power to the external estimate of the mean 
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Figure 3. Gaussian realisations of the ID power spectrum predicted by the 3D power spectrum of the APM survey after convolution 
with the survey window function of Szalay et al. (1993). The vertical dotted line shows the wavenumber of the spike in the power 
spectrum reported by BEKS and Szalay et al. (1993). The lower dashed line shows the mean power level at k = 0.049/iMpc^ 1 . The 
upper dashed line shows the peak amplitude reported by Broadhurst et al. (1995). Two realisations are shown (a) and (c) that contain 
a sharp peak with no other peak higher than P ~ 0.02; such spectra occur much less often than the type shown in (b), where there are 
peaks several times higher than the mean power at other wavenumbcrs. Panel (d) shows and example of a spectrum without any high 
peaks. 



power. For the BEKS data, (geometries A and B), the peak 
power is P = 80.92. KP point out that 380 galaxies have 
redshifts less than z = 0.5, which gives an amplitude in our 
units of P = 80.92/380 = 0.213, which is used to calculate 
the peak to mean power ratio for cases A and B. Figure 4 
of Broadhurst et al. (1995) gives the height of the peak for 
many beams as P = 0.085, which is used in case C. 

KP found that projected power spectra have an expo- 
nential distribtuion of amplitudes, as expected from the cen- 
tral limit theorem (Fan & Bardeen 1995): the probability 
of observing a peak in the power spectrum in excess of P 
at some selected frequency given the mean power (P) is 
(Kendall & Stuart 1977) Pr(> P) = exp(-P/(P}). 

There is a factor of ~ 2 difference in the peak to mean 
power ratio for different survey geometries which illustrates 
the difficulty in making an external estimate of the mean 
power. When the effects of redshift space distortions are in- 



cluded, the peak to mean power ratio can fall as low as 5.7, 
which gives a field point probability that the peak arises 
from a Gaussian random field of 3.3 x 10~ 3 . This is much 
higher than the most pessimistic example in which there are 
no redshift space distortions and a = 0, where the proba- 
bility is 4.5 x 10~ 5 . However, in both cases the field point 
probability is much larger than the figure 4 x 10~ 7 quoted by 
Broadhurst et al. (1995) for the null hypothesis of no large 
scale clustering of galaxies. Note that the amplitude of the 
APM power spectrum could be low by as much as 20%, due 
to merging corrections to APM images which we have not 
applied here (Maddox et al. 1990, 1996) or uncertainties in 
the form of the redshift distribution (Baugh & Efstathiou 
1993, Gaztanaga 1995). 

We have simulated ID power spectra using the pro- 
jected amplitude for Pic(fe) computed from the APM power 
spectrum after convolution with the multiple beam win- 
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dow function. Real and imaginary Fourier components are 
generated with a Gaussian distribution. The resolution of 
the ID power spectrum is taken to be 8k = 0.007hMpc _1 
as calculated by KP for the combined sample of BEKS 
and which matches the full width at half maximum of the 
peak in Figure 4 of Broadhurst et al. (1995). Examples of 
these synthetic spectra are shown in Figure |j| We gener- 
ate a large number of realisations of the Gaussian power 
spectrum and look for peaks in the wavenumber range 
0.028hMpc _1 < k < 0.063/iMpc -1 , corresponding to wave- 
lengths of 100/i,~ 1 Mpc < A < 224/i _1 Mpc. We find peaks 
in excess of the height reported by Broadhurst et al. (1995) 
3% of the time using the estimated real space mean power, 
rising to 6% of realisations when redshift space distortions 
with (3 = 1.0 and a v = 500kms -1 are included. A typical 
power spectrum generated with a high peak in the range 
reported by Broadhurst et al. is shown in Figure ^(b). The 
power spectrum realisation contains several peaks that are 
a few times higher than the mean power, indicated by the 
dotted line, in contrast with the power spectrum shown by 
Figure 4 of Broadhurst et al. . Roughly 1 in every 1500 re- 
alisations contains a peak in excess of P = 0.085 in the 
specified wavenumber range, without having another peak 
in the range 0.0 < k < 0.3 that is higher than P = 0.02; two 
examples are shown in Figure |(a) and (c). 



5 DISCUSSION 

We have examined the compatibility of the three dimen- 
sional power spectrum measured from the APM Galaxy Sur- 
vey (Baugh & Efstathiou 1993, 1994) with the power spec- 
trum recovered from deep pencil beam surveys by BEKS, 
Szalay et al. (1993) and Broadhurst et al. (1995). The prob- 
ability of observing a peak in the ID power spectrum, given 
our knowledge of the 3D power spectrum on these scales 
depends upon the ratio of the height of the observed peak 
compared with the external estimate of the mean power. Ex- 
ternal estimates of the mean power in the ID datasets are 
extremely sensitive to how the survey geometry is modelled 
and to the magnitude of the redshift space distortions. For a 
survey consisting of many pencil beams, we find a factor of 
two difference in the amplitude of the mean power that we 
calculate, for various models of the distortions in the pattern 
of clustering due to the peculiar velocities of clusters. 

In the case of our highest estimate of the mean power 
for a survey consisting of multiple pencil beams, we find that 
the field point probability of finding a peak of the observed 
height given the null hypothesis of the measured clustering 
in 3D, is roughly 1 in 300. We have also generated Gaussian 
realisations of the pencil beam power spectra, using the ID 
power spectrum obtained by convolving the measured 3D 
power spectrum with the survey geometry. Using the same 
resolution as the observed ID power spectrum, we find that 
peaks of the observed height or higher in the wavelength 
range 100 < A < 200/t" 1 Mpc occur in 3 - 6% of the reali- 
sations, depending upon the model assumed for the redshift 
space distortions. 

The general appearence of these synthetic power spectra 
is quite different from that of the power spectrum of the 
observations in the majority of realisations. We find that 
a high peak in the power in the wavelength range specified 



above is generally accompanied by peaks that are a few times 
the amplitude of the mean power at other wavenumbers. In 
only a small number of cases, roughly 1 in 3000 do with find 
an isolated high peak with the other peaks in the spectrum 
less than twice the amplitude of the mean. KP discuss the 
effects that redshift errors or binning of the data in real 
space before taking the Fourier transform could have upon 
the appearence of the observed power spectrum. Also, we 
have neglected the consequences of the density distribution 
in 3D being skewed (e.g. Efstathiou et al. 1990, Saunders et 
al. 1991, Gaztahaga 1994) which has been shown to further 
enhance the probability of observing a peak given a null 
hypothesis of the measured 3D clustering (Amendola 1994). 

The ability of the measured 3D clustering to account 
for the presence of high peaks in the observed ID power 
spectrum makes it seem unlikely that large scale structure 
can be responsible for the steep observed slope of the galaxy 
counts at bright magnitudes, a conclusion already discussed 
by Maddox et al. (1990) and Loveday et al. (1992). The 
factor of two increase in the normalisation of the luminosity 
function desired to provide a better match to the faint counts 
(Shanks 1990) requires that the local universe out to a radius 
of around 130fe _1 Mpc, corresponding to the median redshift 
of galaxies brighter than bj = 17, is massively underdense. 
The variance in the number of galaxies measured on these 
scales from the APM Galaxy Survey is on the order of 10% 
(Baugh & Efstathiou 1993, Gaztanaga 1994) many times 
smaller than would be necessary to result in a factor of two 
fluctuation in the number of galaxies found in a volume of 
such a large radius (see also the discussion in Glazebrook et 
al. 1994). 
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